arXiv:1502.01030v2 [hep-ph] 30 Apr 2015 


DESY-15-013, TTP15-01, SFB/CPP-14-120 


Quark mass relations to four-loop order in perturbative QCD 

Peter Marquard, 1 Alexander V. Smirnov, 2 Vladimir A. Smirnov, 3 and Matthias Steinhauser 4 
1 Deutsches Elektronen Synchrotron DESY, Platanenallee 6, 15738 Zeuthen, Germany 

2 Scientific Research Computing Center, Moscow State University, 119991, Moscow, Russia 

3 Skobeltsyn Institute of Nuclear Physics, Moscow State University, 119991, Moscow, Russia 
4 Institut fur Theoretische Teilchenphysik, Karlsruhe Institute of Technology (KIT), 76128 Karlsruhe, Germany 

(Dated: February 2015) 

We present results for the relation between a heavy quark mass defined in the on-shell and MS 
scheme to four-loop order. The method to compute the four-loop on-shell integral is briefly described 
and the new results are used to establish relations between various short-distance masses and the MS 
quark mass to next-to-next-to-next-to-leading order accuracy. These relations play an important 
role in the accurate determination of the MS heavy quark masses. 

PACS numbers: 12.38.Bx, 12.38.Cy, 14.65.Fy, 14.65.Ha 


The precise knowledge of quark masses plays an im¬ 
portant role in many phenomenological applications. 
This is in particular true for the heavy top, bottom and 
charm quarks. For example, the top quark mass enters as 
a crucial parameter the combined electroweak fits which 
have been used to obtain indirect information about the 
Higgs boson mass, and nowadays serve as consistency 
checks for the Standard Model, see, e.g. Refs. M- The 
uncertainty in the top quark mass is also dominant in the 
analyses of the stability of the electroweak vacuum 0- 
E]. A prominent example where a precise bottom quark 
mass is required are R-meson decays which are often pro¬ 
portional to the fifth power of mb- Precise charm and 
bottom quark masses are important to obtain accurate 
predictions for the Higgs boson decays into the respective 
quark flavours. Also in the context of top and bottom 
Yukawa coupling unification, precise mass values are in¬ 
dispensable since they serve as boundary conditions at 
low energies. Last but not least, quark masses enter the 
Lagrange density of the Standard Model as fundamental 
parameters. Thus, it is mandatory to obtain precise nu¬ 
merical values by comparing high-order theoretical pre¬ 
dictions with precise experimental data. 

At lowest order in perturbation theory there is no need 
to fix the renormalization scheme for the quark masses. 
However, after including quantum corrections it is neces¬ 
sary to apply renormalization conditions which fix the 
renormalization scheme. A natural scheme for heavy 
quark masses, i.e. the charm, bottom and top quark 
masses, is the on-shell (OS) scheme where one requires 
that the inverse heavy quark propagator with momentum 
q has a zero at the position of the on-shell mass, M, i.e. 
for q 2 = M 2 . It is well known that perturbation the¬ 
ory has a bad convergence behaviour in case the on-shell 
quark mass is used as a parameter. Another widely used 
renormalization scheme is based on minimal subtraction. 
This means that the mass parameter entering the quark 
propagator is defined in such a way that just divergent 
terms (and no finite contributions) are absorbed such 
that the quark propagator is finite (after wave function 
renormalization). In this Letter we consider four-loop 
corrections to the relation between the on-shell and the 


MS definition of a heavy quark mass which allows for a 
precise conversion from one renormalization scheme into 
the other. 

For the various heavy quarks different methods relying 
on different quark mass definitions are used to extract the 
mass values. For example, in Ref. [f| low-moment sum 
rules have been used to extract directly the MS charm 
and bottom quark masses without any reference to the 
on-shell mass. On the other hand, physical observables 
inherently connected to the threshold, like T sum rules 
or top quark pair production close to threshold, rely on 
properly defined quark masses, like the potential subtrac¬ 
ted (PS) 0, IS (8 hIo| or renormalon subtracted (RS) [ll| 
definition. When comparing with experimental data, in 
a first step the corresponding mass values are extracted. 
Afterwards they are converted to the MS definition. Note 
that the relation between the MS and the OS mass is an 
important ingredient to obtain the relation between the 
PS, IS or RS masses and the MS mass. 

In this Letter we use the four-loop MS-OS relation to 
establish relations between the PS, IS, RS and the MS 
quark mass which are necessary to obtain the latter with 
next-to-next-to-next-to leading order (N 3 LO) accuracy. 

Note that there is a further definition of a threshold 
mass, the so-called kinetic mass [l2| which has been used 
for quite a number of applications in B physics, (see, e.g., 
Ref. M)- However, the relation to the on-shell mass is 
only known to two loops (NNLO). For this reason it is 
not considered in the following. 

In the following we first discuss the relation between 
the MS and OS quark mass. Afterwards we elaborate on 
the relation between the threshold (PS, IS and RS) and 
the MS mass. The latter is obtained by using as starting 
point the definition of the PS, IS or RS masses which 
establishes a relation to the pole mass. Afterwards the 
pole mass is replaced by the MS mass which leads to the 
desired relation between the short-distance masses. 

To derive a formula relating the MS and OS quark 
mass it is advantageous to start with relations between 
theses masses and the bare mass, to 0 , which are given 
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b>0 

m° = Z^m , m° = Z° S M . (1) 

Z^ s is known to four loops and can be found in Refs. [fi¬ 
lial . By construction, the ratio of the two equations in © 
is finite which leads to 

2m (w = > (2) 

where z m depends on a s (/j,) and log (fi/M) and has the 
following perturbative expansion 

= ( 3 ) 

n> 0 


with Zm = 1. 

For a derivation of convenient formulae relating 
to the on-shell quark self energy we refer to Refs, 
where it is shown that Z^ s is obtained from the sum of 
scalar and vector contribution evaluated on-shcll, i.e., 

Z° S = 1 + S y (g 2 = M 2 ) + £ s (g 2 = M 2 ). (4) 



where f3o = 11 — 2n;/3 is the one-loop coefficient of the 
QCD /? function and a± = 31/3 — 10rq/9 the one-loop 
coefficient of the static potential, ni is the number of 
massless quarks, is the factorization scale which is of 
the order of the soft scale. In this Letter we use f-if = 
2 GeV for bottom and fif = 20 GeV for top quarks. 
Sm(fif) is known to N 3 LO ||28] which involves the three- 
loop corrections to the static potential, 03 [29M3TI ] . 

The N 3 LO relation between ro PS and m is obtained by 
inserting Eq. (j5j into Eq. (j6|). All ingredients are already 
expanded and the coefficients of (a s ) n have simply to be 
combined. In particular, the term in ([5]) involving 03 is 
combined with the four-loop term in the m-M relation. 
One obtains an explicit formula to compute the PS mass 
in case the MS is given. For a given PS mass we solve 
this equation iteratively to obtain the MS mass. 

The IS mass is defined as half the perturbative mass of 
a fictitious l 3 Si state, where it is assumed that the quark 
is stable. Thus, we have the following relation between 
the IS and on-shell mass [m3 


TO 1S = M + -E 




(9) 


One-, two- and three-loop QCD results to Z° s have been 
computed in Refs. 0], [l7f and [Hi EH HH S3 , respect¬ 
ively, and electroweak effects have been considered in 
Refs. [H - l27j . The main task of this Letter is the com¬ 
putation of the four-loop QCD corrections to Z^ s and 
consequently to z m . For convenience we introduce also 
the inverse relation to Eq. m as follows 

M = j). (5) 

The PS quark mass has been introduced in Ref. (Tf. 
Its relation to the pole mass is given by 

m FS = M — 5m(iJ,f), (6) 


with 


Sm(iif) 


1 

2 


'l<? I <Hf 


d 3 g 

(2/Q 3 


V(q), 


(7) 


where V(q) is the perturbative contribution to the static 
heavy quark potential. can be computed in per¬ 

turbative QCD and has the form 





( 8 ) 


1 We refrain from adding a superscript “OS” to the on-shell mass 
but use a capital letter. Similarly, a lower-case m without further 
superscript stands for the MS quark mass. For all other mass 
definitions we use a lower-case “m” and a superscript indicating 
the renormalization scheme. 


where is the perturbati ve g round state energy which 
is available to third order |2T [32| , [33| . The last miss¬ 
ing ingredient was the three-loop static potential which 
has been evaluated in Refs. [29l - [3ll ] . The replacement 
a" —> a"e n_1 implements the so-called e expansion 
which guarantees that the appropriate orders in the ex¬ 
pansions of and M (in terms of m) are combined. 

The perturbative expansion of E pt has the following 
form 

Ef 

“ " n> 0 


UhMaz 


= -£- 


E 


(s^ 
\ 7r 


s:( n ) 


( 10 ) 


In order to obtain the relation between m is and m one 
has to replace in Eq. fl5J by a^ ll ' > and then ap¬ 

ply the replacement a" —>• a"e n . Afterwards it is inser¬ 
ted into Eq. (|!J]) and expanded in the parameter e. This 
guarantees that the N fe LO term in is combined with 
the (k + l)-loop correction to c m (/x). In particular, in 
order to establish the m is -m relation to N 3 LO four-loop 
corrections to c m (/i) are needed. A given MS quark mass 
is transformed to the IS mass by inserting the numerical 
value into the resulting equation. In case the IS mass is 
given one obtains the MS mass by solving the equation 
implicitly. 

A further threshold mass, the so-called RS mass, has 
been introduced in Ref. & It is related to the pole mass 
in such a way that the pure renormalon contributions 
are subtracted. The corresponding formulae are derived 
and explicitly given in Ref. m- In that reference also 
a variant, the so-called RS' scheme, is discussed where 
no subtraction is performed for the 0(a s ) term in the 
MS-OS relation. Recently the numerical accuracy of the 
normalization constant of the first renormalon has been 
improved in Ref. [3dj , where also a variant of the RS and 
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RS' masses has been suggested in which (in the case of the 
bottom quark) the subtraction term is parameterized in 

(3) 

terms of a), . In this paper we will adopt the prescription 
of Ref. Q Similarly to the PS mass also for the RS mass 
a subtraction scale has to be specified which we again 
choose &s fif = 2 GeV for bottom and pf = 20 GeV for 
top. 

For the computation of the scalar and vector part of 
the fermion propagator we use an automated setup which 
generates all contributing amplitudes, processes them 
with F0RM3 0] and provides scalar functions involving 
several million different integrals encoded in functions 
with 14 different indices which belong to 100 different 
integral families. 

The Laporta algorithm 0] is applied to each family 
using FIRE5 [37 ] and crusher l3jj w hich are written in 
C++. Then we use the code tsort [39|, which is part of the 
latest FIRE version, to reveal relations between primary 
master integrals following recipes of [40 ] and end up with 
386 four-loop massive on-shcll propagator integrals, i.e. 
with p 2 = M 2 . 

We have performed the calculation allowing for a gen¬ 
eral gauge parameter £ keeping terms up to order £ 2 in 
the expression we give to the reduction routines. We have 
checked that £ drops out after mass renormalization but 
before inserting the master integrals. 

For some master integrals, analytic results could be de¬ 
rived using a straightforward loop-by-loop integration for 
general space-time dimension. We also used analytical 
results obtained for non-trivial four-loop on-shcll mas¬ 
ter integrals computed in our earlier paper Ref. [4l| . In 
some other cases one- and two-fold Mellin-Barnes repres¬ 
entations can be derived which allow for a high-precision 
numeric evaluation, at least up to 20 digits. For some 
of the master integrals, we applied threefold MB repres¬ 
entations which enabled us to obtain a precision of eight 
digits. 

For factorizable integrals, we obtained analytic results 
from known two- and three-loop results. In particular, 
we used Ref. 0 ] where the expansion in e = (4 — d)/2 
has been performed up to the order typical to four-loop 
calculations. (d is the space-time dimension used to com¬ 
pute the momentum integrals.) 

We computed the remaining 332 integrals numerically 
with the help of FIESTA [43l - l45| . FIESTA returns for each 
e coefficient a numerical result and the corresponding un¬ 
certainty from the numerical integration. When inserting 
the master integrals we keep track of all uncertainties and 
combine them quadratically in the final expression. We 
interpret the resulting uncertainty as a standard devi¬ 
ation and multiply it by five in the final result for the 
relation between the MS and OS quark mass. This is in 
agreement with adding the uncertainties from the indi¬ 
vidual contributions linearly. 

We are now in the position to present numerical res¬ 
ults for z m {p) which have been obtained by setting the 
number of colours to three (N c = 3) and the number of 
massless quarks ( ni ) to either 3, 4 or 5, corresponding to 


the charm, bottom or top quark case, before combining 
the uncertainties from the numerical integration of the 
master integrals. Note that the coefficients up to three 
loops are known analytically 0 , 0 ]. We refrain from 
listing the corresponding results but refer to Eq. (13) of 
Ref. 0] . Analytical results are also available for the log¬ 
arithmic four-loop contributions since they can easily be 
obtained using renormalization group methods. In the 
following we restrict ourselves to compact numerical res¬ 
ults. At four loops we obtain for the coefficient of (a s /ir) A 


z 


(4) 


ni— 3 


Z 


(4) 


ni=4 


-1744.8 ±21.5 - 703.48 Z OS 
- 14.234 Z^ s - 0.75043 Z£ s , 
-1267.0 ± 21.5 - 500.23 Z os 


122.97 Iq S 

83.390 Zq S 


z 


(4) 


ni— 5 


- 9.9563 Z^ s -0.514033 Z^s 
-859.96 ± 21.5 - 328.94 Zqs 


50.856 Zq S 


- 6.4922 Z^ s - 0.33203 Z^ s , (11) 

with Zos = In {p 2 /M 2 ). We obtain the ^.-independent 
coefficients with an accuracy of 1.2% for rq = 3, 1.7% for 
ni = 4) and 2.5% for m = 5. In the numerical results 
discussed below we will assume a relative uncertainty of 
3% for all values of m. 

For convenience we also show the four-loop results for 
c m which read 


c 


(4) 


ni— 3 


C 


(4) 


rt;=4 


C 


(4) 


ni=5 


1691.2 ± 21.5 + 828.43 l m + 189.65 Z^ 


+ 36.688 Z|js + 4.8124 Z^, 


MS ’ 


1224.0 ± 21.5 + 601.98 l m + 134.10 Z^ 
+ 28.846 1^+ 3.9648 Z^, 

827.37 ± 21.5 + 408.88 + 86.574 Z^ 


+ 22.023 Z^ + 3.2227 Z 4 (12) 


with Zjyfg = In (p 2 /m 2 ). In the remaining part of this 
Letter we will concentrate on the top and bottom quark 
mass. 

As an application of the new results in Eqs. ED 
and ED we study the relations between the various 
threshold masses and the MS mass. We use the follow¬ 
ing input values for the strong coupling constant and the 
bottom and top quark masses: 

af\M z ) = 0.1185 [13 , m b {m b ) = 4.163 GeV @ , 
M t = 173.34 GeV 0 . (13) 


(5) 

a s with four and six active flavours is obtained from as 
where for the decoupling scale we choose twice the heavy 
quark mass 0,0 . 

Let us have a closer look to the relation between the 
OS and MS top quark mass. For p = m t we have 


M t = m t (1 + 0.4244 a s + 0.8345 a 2 s + 2.375 a 3 s 
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+(8.49 ± 0.25) a*) 

= 163.643 + 7.557 + 1.617 + 0.501 

+ 0.195 + 0.005 GeV, (14) 

with a s = ai 6 \m t ) = 0.1088. Note that the four-loop 
term still gives a contribution of about 200 MeV which 
is not negligible even with nowadays uncertainties from 
TEVATRON and LHC [5j§]. The corresponding results 
for the bottom quark read 

M b = m b ( 1 + 0.4244 a s + 0.9401 a 2 s + 3.045 a 3 s 
+(12.57 + 0.38) afj 
= 4.163 + 0.401 + 0.201 + 0.148 

+ 0.138 + 0.004 GeV. (15) 

Here a s = ai 5 \m b ) = 0.2268. Note that the four-loop 
corrections in Eq. © are almost as large as the three- 
loop term. On the other hand, the perturbative series 
for the case of the top quark has a reasonable behaviour: 
the three-loop coefficient is by a factor three smaller than 
the two-loop one and the four-loop term is again smal¬ 
ler by a factor 2.5. This suggests that with the help of 
Eq. (1141) the top quark mass can be determined with an 
uncertainty below 200 MeV. 

In practice it often happens that in a first step a 
threshold quark mass is extracted from comparisons 
of higher order calculations and experimental measure¬ 
ments. Afterwards the threshold mass is converted to the 
MS quark mass. In Tabs. [Q and El we show the results 
for the scale invariant MS quark mass m q (rn q ) ( q = b, t) 
using one- to four-loop accuracy for the conversion. 


input 

m PS = 

m ls = 

m RS = 

Sloops 

171.792 

172.227 

171.215 

1 

165.097 

165.045 

164.847 

2 

163.943 

163.861 

163.853 

3 

163.687 

163.651 

163.663 

4 

163.643 

163.643 

163.643 

4 (x 1.03) 

163.637 

163.637 

163.637 


Table I. mt(mt) in GeV computed from the PS, IS and RS 
quark mass using one- to four-loop accuracy. The numbers 
in the last line are obtained by increasing the four-loop coef¬ 
ficient in Eq. (THU) by 3%. 

In the case of the top quark (cf. Tab. Q} the three- 
loop corrections amount to about 200-250 MeV which re¬ 
duces to {44,8,20} MeV at four loops for the {PS,1S,RS} 
quark mass. A 3% uncertainty in the MS-OS relation 
induces a shift of 6 MeV in mt(mt) which is in gen¬ 
eral small as compared to the four-loop contribution. 
Let us estimate the final uncertainty from the conver¬ 
sion to the MS mass from the quadratic combination 
of the 6 MeV with half of the four-loop contribution 
(i.e. (44,8 ,20}x 1/2 MeV). This leads to {23,7,11} MeV 
which should be added in quadrature to the remaining 
uncertainties of the threshold mass. 


input 

#loops 

m PS = 
4.483 

m ls = 

4.670 

m RS = 

4.365 

1 

4.266 

4.308 

4.210 

2 

4.191 

4.190 

4.172 

3 

4.161 

4.154 

4.158 

4 

4.163 

4.163 

4.163 

4 (x 1.03) 

4.159 

4.159 

4.159 


Table II. m b (mb) in GeV computed from the PS, IS and RS 
quark mass using one- to four-loop accuracy. The numbers 
in the last line are obtained by increasing the four-loop coef¬ 
ficient in Eq. (1121) by 3%. 


The results for m b (m b ) computed from the PS, IS and 
RS threshold masses are shown in Tab. El The three-loop 
corrections provide still sizable effects of up to 40 MeV 
which reduces to at most 9 MeV at four loops. The uncer¬ 
tainty in the four-loop MS-OS relation induces an error 
of 4 MeV. Thus we arrive at a final error of {4, 6, 5} MeV 
for the conversion from the {1S,PS,RS} mass. This is 
not negligible, though in general much smaller than other 
uncertainties involved in the quark mass extraction (see, 
e.g., Refs. d3|, Q an d [U] for recent determinations of 
TO ft( TO &) where in intermediate steps the IS, RS and PS 
has been used, respectively). 

The results of Tabs. Q] and El can be used, in com¬ 
bination with similar calculations for different values of 
a s (Mz) and threshold masses, to construct the following 
approximation formulae 

= 163.643 ± 0.023 + 0.074A Qs - 0.095A PS , 
GeV t 

= 163.643 ± 0.007 + 0.069A™ - 0.096A1® , 
GeV 3 m * 

= 163.643 ± 0.011 + 0.067A a - 0.095A?. S , 
GeV * 

= 4.163 ± 0.004 + 0.007A a - 0.018A?® , 

GeV b 

= 4.163 ± 0.006 + 0.008A™ - 0.019A^ , 

GeV b 

= 4.163 ± 0.005 + 0.004A as - 0.018A?f (16) 

GeV 

with A Qs = (0.1185 - a,(M z ))/0.001, A ps ( = 

(171.792 GeV - mf s )/0.1, A)„ s t = (172.227 GeV - 
m t is )/0.1, A™ = (171.215 GeV - mf s )/0.1, A PS b = 
(4.483 GeV—m PS )/0.02, A^ 6 = (4.670 GeV-mf )/0.02, 
A* s b = (4.365 GeV - mf s )/0.02. 

Let us finally compare in Tab. Mil our result for 
the four-loop coefficient c>m to predictions obtained on 
the basis of different assumptions. In general good 
agreement is found, in particular with the results from 
Refs. [34l . f55l . [56 ] which are all based on renormalon can¬ 
cellation. For example, in Ref. [56j . the four-loop coef¬ 
ficient is extracted from the requirement of perturbative 
stability of the combination 2m po i e + Vqcd where Vqcd 
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Ref. 

= m(m)) 


ni = 3 

ni = 4 

ni = 5 

[52] 

[53] 

M 

[55] 

[56] 

m 

1668 

1571.4 

1281 

1785.9 

1668 ± 167 

1772 ± 82 

1324 

1107.8 

986 

1316.4 

12581” 

1324 ± 82 

1031 

727.0 

719 

920.1 

897t?75 

945 ± 92 

this work 

1691.2 ±21.5 

1224.0 ±21.5 

827.37 ±21.5 


Table III. Comparison of predictions for the four-loop coeffi¬ 
cient Cm' 1 11 to our result which is present in the last line. 


is the static potential of two heavy quarks. The estimates 
in Ref. [53] have been obtained only on the basis of the 
two-loop results, leading nevertheless to good approxima¬ 
tions. Somewhat lower results have been obtained in 
Ref. [54j] where dispersive methods have been used and 
large n 2 terms have been identified. In Ref. [12 the four- 
loop relation between the on-shell and MS quark mass 
has been estimated using the large /3q approximation. 

To conclude, in this Letter we have computed the four- 


loop corrections between the on-shell and MS definition 
of heavy quarks. Our main results are given in Eqs. m 
and (1121) for charm, bottom and top quarks. As applic¬ 
ations we have derived precise relations between the PS, 
IS and RS threshold masses and the MS quark mass. 
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